Decay of superfluid currents in the interacting one-dimensional Bose gas 
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We examine the superfluid properties of a one-dimensional (ID) Bose gas in a ring trap based on the model 
of Lieb and Liniger. While the ID Bose gas has nonclassical rotational inertia and exhibits quantization of ve- 
locities, the metastability of currents depends sensitively on the strength of interactions in the gas: the stronger 
the interactions, the faster the current decays. It is shown that the Landau critical velocity is zero in the thermo- 
dynamic limit due to the first supercurrent state, which has zero energy and finite probability of excitation. We 
calculate the energy dissipation rate of ring currents in the presence of weak defects, which should be observable 
on experimental time scales. 
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I. INTRODUCTION 

Superfluidity is one of the most dramatic manifestations of 
quantum mechanics on the macroscopic scale, and is associ- 
ated to a host of different phenomena such as non-classical 
rotational inertia, quantization of vortices, dragless motion of 
impurities and metastability of ring currents as seen in, e.g., 
liquid He II. Since each of these phenomena may be taken 
as "defining" a transition to superfluidity, it is important to 
ask under what circumstances they occur together. As was 
pointed out by Leggett lUt] the metastability of ring currents 
and nonclassical rotational inertia are two fundamental super- 
fluid phenomena of yet very different nature. While the latter 
is an equilibrium property, the former is a dynamic one. Al- 
though both types of phenomena are often explained by Bose- 
Einstein condensation of bosons or Cooper pairs of fermions 
lUt], the latter is not seen as an exclusive requirement 10, Ht]. 
Here we consider the superfluid properties of an interacting 
one-dimensional (ID) Bose gas at zero temperature, a system 
which is not Bose-condensed |0, Ht] but may possess quasi- 
long -range order |6fl. It is a long-standing question whether 
the ID Bose gas can support persistent currents with macro- 
scopic lifetimes jst]. 

This system has been realized with ultracold bosonic atoms 
in tightly confining linear traps 10, 01 (ring traps are also un- 
der development ||9|])jin which the boson interactions are ef- 
fectively described ITol [Till by the contact potential V{x) = 
gBS{x) of the Lieb-Liniger (LL) model lll2ll . The interac- 
tion strength is quantified by the dimensionless parameter 
7 = mg-Q/ {h^n), where n is the linear density and m is 
the mass. For 7 00, the model is known as the Tonks- 
Girardeau (TG) gas and can be mapped to an ideal Fermi 
gas. For 7^1, the Bogoliubov model of weakly interact- 
ing bosons is recovered. 

Experimental investigation of the superfluid properties of 
the ID Bose gas by observing the motion of impurities is at an 
early stage and theoretical predictions are not yet compre- 
hensive. Sonin |0] found that ring currents can be metastable 
except for infinitely strong interactions. Kagan et at. lfl3ll also 
concluded that persistent currents could be observable on ex- 



perimental time scales and Biichler et al. ifTill found the ID 
Bose gas able to sustain supercurrents even in t he p resence 
of a strong defect. Astrakharchik and Pitaevskii lUSh consid- 
ered the drag force on a moving heavy impurity within Lut- 
tinger liquid theory and predicted a power-law dependence on 
the velocity for small velocities. These results contain an un- 
known prefactor preventing the calculation of the actual value 
of the drag force and are in any case not applicable at larger 
velocities. The motion of an impurity of finite mass was con- 
sidered in the TG gas [16] but for finite values of 7 this prob- 
lem is still unresolved. 

In this paper we calculate the rate of energy dissipation of 
ring cuiTents in the presence of a small integrability-breaking 
defect of strength q,. based on recent advances in the under- 
standing IitI [isl [1911 of the dynamics of the LL model. The 
results of our calculations are summarized in Fig. [T] While 
for small velocities our calculations support the power-law 
predictions of Ref. ifTsll . the drag force = 2gfnmhr^ 
assumes the velocity-independent value of 2gfnm/h^ for ve- 
locities large compared to the speed of sound c. Although our 
results suggest that the ID Bose gas can support metastable 
currents only in the weakly interacting regime where 7^1, 
the superfluid fraction is 1 regardless of 7 ll20ll according to 
the nonclassical rotational inertia for a finite ring. 

II. LANDAU CRITERION OF SUPERFLUIDITY 

In the LL model the total momentum is a good quantum 
number, and periodic boundary conditions quantize it in units 
of 2TTh/L, where L is the ring circumference. The low-lying 
spectrum of iV = nL bosons as shown in Fig.|2]has local min- 
ima ll2lll at the supercurrent states / (/ = 0, 1, 2, . . .) with mo- 
menta p/ = 2'iTnhI and excitation energies £/ = pj/ {2Nm). 
These correspond to Galilean transformations of the ground 
state with velocities vj = pj/{Nm). The minima do not 
depend on interactions and tend to zero in the limit of large 
system size. 

Suppose that the gas is initially rotating with the linear ve- 
locity —vj and then is braked with an "obstacle," created, e.g.. 
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FIG. 1: (Color online) The dimensionless drag force versus the ve- 
locity (relative to vf ~ h-Kn/m) of the impurity at various values 
of the coupling parameter. The solid (blue) lines represent the force 
obtained with Eqs. l|T) andjRJ; open circles are the numerical data 
obtained using ABACUS (m . 




FIG. 2: (Color online) Schematic of the excitation spectrum of the 
ID Bose gas in a perfectly isotropic ring. The supercurrent states I 
lie on the parabola h^k^ / {2Nm) (dotted line). Excitations occur in 
the shaded area; the discrete structure of the spectrum is not shown 
for simplicity. The blue (dark) area represents particle-hole excita- 
tions 12311 . Motion of the impurity with respect to the gas causes 
transitions from the ground state to the states lying on the straight 
(red) line. 

by a laser beam ll22ll . In the frame where the gas is at rest, the 
obstacle moves with velocity vi. In a superfluid we expect 
no energy dissipation, and thus zero drag force (the current is 
persistent). Energy conservation dictates that the transitions 
from the ground state caused by the moving obstacle with ve- 
locity V, lie on the line e = vp. According to Landau, if 
the excitation spectrum lies above this line, the motion can- 
not excite the system, which is then regarded superfluid. The 



Landau critical velocity (when the line touches the spectrum) 
equals = ei/pi ~ fi/2. This implies that any supercurrent 
state with / ^ 1 is unstable since vj > Vc- However, in 3D 
similar supercurrent states exist, which apparently leads to the 
absence of current metastability. The paradox can be resolved 
by considering not only the spectrum but also probabilities of 
excitations. Below we argue that in the 3D case, the probabil- 
ity to excite supercurrents is vanishingly small, while in the 
ID case it depends on the strength of bosonic interactions. 

III. HESS-FAIRBANK EFFECT 

When the walls of a toroidal container are set in rotation 
adiabatically with a small velocity, a superfluid stays at rest 
while a normal fluid follows the container This effect leads to 
a nonclassical rotational inertia of superfluid systems, which 
can be used to determine the superfluid fraction One can 
show ll20ll that in the LL model, the gas has zero rotational 
inertia (zero normal fraction) at T = for any 7 > 0. This 
is an equilibrium property completely determined by the low- 
lying energy spectrum ji^ . 

IV. DECAY OF SUPERCURRENTS 

A. Dynamic response 

By contrast to the Hess-Fairbank effect, metastability of 
currents is not an equilibrium effect and transition probabil- 
ities have to be considered. The dissipation rate as energy loss 
per unit time E of an obstacle (or heavy impurity) moving 
with velocity v relative to the gas can be related to the drag 
force Fv acting on the impurity hy E = —F^v. For weak im- 
purities with interaction potential Vi{x) the drag force is re- 
lated to the dynamic structure factor (DSF) in linear-response 
theory d Hi: 

F^{v)= dkk\Vi{k)\^S{k,kv)/L, (1) 

Jo 

where Vi (fc) is the Fourier transform of the impurity potential. 
The DSF S{k, uj) describes the transition probability between 
the ground state |0) and excited states \ m) with energy transfer 
huj and momentum transfer Hk caused by a density perturba- 
tion, and can be written as 

S{k,Lo) = ^ |(0|,5p;,|m)|2^(;iw - E„, + Eo), (2) 

m 

where Spk = J2j e"*'^^^ — NA{k) is the Fourier component 
of the density operator, A(fc) = 1 at /c = and A{k) = 
otherwise. Several results for the DSF in the LL model have 
recently become available lfl7l[l8l[l9t| . It can be measured in 
cold gases by Bragg scattering ll25ll26ll . 

Numerical values of DSF calculated with the ABACUS al- 
gorithm ifisll are shown in Fig. [3] The probability to create 
multiparticle excitations lying outside of the region (fc) ^ 
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FIG. 3: (Color online) Dynamic structure factor of the ID Bose gas 
from llal for A'^ — 100. Dimensionless values of S{k,uj)£F/N 
are shown in shades of gray between (white) and 0.7 (black). The 
full (blue) lines represent the limiting dispersion relations uj± and 
the straight (red) line is the line of integration in Eq. l[T). Only one 
point at fc = fee, shown in full (red) circle, contributes to the integral 
when the perturber is a shallow cosine potential with a reciprocal 
vector kc- 

u! ^ are identically zero (below ti>_) or very small 

(above oj+). Transitions from the ground state caused by a 
moving obstacle with velocity v occur along the straight (red) 
line. Drag force ^ is thus a generalization of the Landau 
criterion for superfluidity. Indeed, if the excitation spectrum 
of a generic system Ues above the line uj ~ vk then it is su- 
perfluid; in this case the drag force ([T]i equals zero. The drag 
force thus proves to be fundamental and can be considered as 
a quantitative measure of superfluidity. 



small 7, its non-zero values practically localize in vicinity of 
V = u)+{ka) /kc- As there is no sharp transition from super- 
fluid to isolated phase in ID ll28ll . we can put the threshold 
equal to, say, 0.1 of the characteristic value TrglkcN / {%e-p) 
of the drag force (|3]l. Then we get a phase diagram in the w-7 
plane iH] similar to that of Polkovnikov et al. @|. Note that 
in the latter paper, the superfluidity was examined in terms 
of quantum phase slips IBOll . So, the both quasiparticle and 
quantum phase slip description lead to the same results. 

C. Approximate expression for drag force 

In order to study metastability of the Ith supercurrent state, 
we need to calculate the drag force on an obstacle moving with 
the velocity vj relative to the gas. For large system size the 
supercurrent-state velocities are dense and in the thermody- 
namic limit {N — > 00, n = const) we may consider arbitrary 
velocities. We consider the drag force and decay of currents 
in various regimes. 

We calculate the drag force from Eq. ^ by using the inter- 
polating expression 

S{k,uj) = C{uj°' -co^)"- /{uj^-Lu°')''+ (4) 

for cj_(fc) ^ cj ^ u!+{k), and S{k,uj) — otherwise 
lIsTll . Here, K = fmn/{mc) is the Haldaiie parameter ifTsll . 

and fJ,-{k) are the exact exponents 111911 at the borders 
of the spectrum uj+{k) and tj_(fc), and a = 1 + 1/VK. 
The values of LLj±{k) and ii±{k) can be calculated from the 
coupling constant 7 numerically by the methods outlined in 
Refs. IQ^IHtl- The normalization constant C depends on mo- 
mentum but not on frequency and is determined from the /- 
sum rule duo ijjS{q,ijj) = Nq^/{2m). The expression 
(HJi is applicable for all ranges of the parameters fc, u), and 7 
with increasing accuracy at large 7. A more detailed discus- 
sion can be found in Ref. Isill . 



B. Shallow optical lattices 

Equation ([T]l can be verified experimentally for different 
types of obstacles: for Vi{x) = gi5{x) all the points at the line 
contribute to the drag force, while for Vi (x) = gi^ cos(27ra; / a) 
only one point {kc, kcv) in the k-u) plane does, where feg = 
27r/a is the reciprocal lattice vector (see Fig.O. Indeed, sub- 
stituting the Fourier transform into Eq. ^ yields 

F^ = TiglkGS{kG,kGv)/2. (3) 

The filling factor of the lattice is given by 21111 /kc- Equation 
(|3]l can be exploited even in the case of a cigar-shaped quasi- 
ID gas of bosons at large number of particles, because the 
boundary conditions do not play a role in the thermodynamic 
limit. It gives us the momentum transfer per unit time from 
a mo ving shallow lattice, which can be measured experimen- 
tally luTll . At kc = 2nn, corresponding to the Mott insulator 
state in a deep lattice, and at 7 1, the drag force takes 
non-zero values for arbitrary v ^ uj+{kG) /kc- However, at 



D. Numerical results 

We further restrict ourselves to a (5-function impurity in- 
teraction with Vi{k) ~ gi. Results of integrating Eq. ^ are 
shown in Fig. For large velocities the drag force reaches 
the velocity-independent value of 2gfnm/h^. A characteris- 
tic velocity scale is the speed of sound c, which determines 
the transition from a power-law increase to the velocity in- 
dependent regime. The speed of sound c of the LL model is 
proportional to 7 for small 7 but saturates to the value vp for 
large 7 llT2ll . The numerical DSF as per Ref. ifisll was obtained 
for N = 150 particles (7 = 20, 5), N = 200 (7 = 1) and 
N — 300 (7 = 0.25). The /-sum rule saturations at fc = 2kp 
were 99.64%(7 = 20), 97.81%(7 = 5), 99.06%(7 = 1), and 
99.08%(7 = 0.25), with yet better results at smaller momen- 
tum. The fit with the analytical ansatz is good for all values of 
V for large 7. The decreasing curves at large velocity v ^ c 
are due to imperfect sum-rule saturation at high momenta. For 
small 7, the onset of the drag force is quicker from the numer- 
ical DSF than from the analytical ansatz. This occurs first 
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because the smoothing of the numerical data required to com- 
pute the drag force overestimates it when its curvature is pos- 
itive (this smoothing also leads to small artifacts in the data 
around v = c), and second because the obtained numerical 
DSF is larger than the analytical ansatz for uj <C 0-!+, and 
also just above the Bogoliubov dispersion (where the analyti- 
cal ansatz is zero by definition), where excitations with higher 
numbers of particle-hole pairs contribute. 

E. Drag force at small velocities 

For the important question whether persistent currents may 
exist at all, the small velocity regime is most relevant, which 
is dominated by transitions near the first supercurrent state 
(umklapp point at fc = 2fcF). The drag force in this regime 
has a power-law dependence on the velocity ~ v'^^^^ for 
V c, as first found by Astrakharchik and Pitaevskii. From 
Eqs. ([T]i and (|4|i we can obtain: 

r(l + ^ - AH(2A:f)) r(l + /i_(2A:F))r(l + 

r(^)r(i-AH(2fcF)) r(i + //_(2fcF) + ^) 

(5) 

where r{x) is Euler's gamma-function, and ^-{2kp) = 
2\/K{\/K — 1) ifigll . This formula is valid for arbitrary cou- 
pling constant and works even in the Bogoliubov regime at 
7 <C 1. In practice, Eq. Q works well up to v < 0.1c. 

F. Why excitations near the umklapp point do not play a role 
in three dimension 

The behavior of the DSF near the umklapp point means 
that the drag force takes non-zero values even for arbitrar- 
ily small interactions. This fact is related to the absence of 
Bose-Einstein condensation in the ID Bose gas. For large in- 
teractions, umklapp excitations become readily available and 
provide an avenue for the rapid decay of supercurrents. Lan- 
dau reasoned that this is very implausible in 3D, since it in- 
volves the macroscopic motion of the system, and hence a 
macroscopic number of quasiparticle excitations. Strong cor- 
relations in ID, however, make umklapp excitations easily ac- 
cessible, since they involve only a single fermionic-like quasi- 
particle lfl2ll . 

G. Currents in a ring 

In the presence of an obstacle a ring current can decay into 
supercurrent states with smaller momentum. Starting in one 
of the local minima of the spectrum in Fig. [2] the kinetic en- 
ergy of the center of mass will be transformed into elemen- 
tary excitations above a lower supercurrent state, conserving 




FIG. 4: (Color online) Decay of the ring current velocity of ID 
bosons from the initial velocity of l.lijp at t = 0. The solid (blue) 
and dashed (red) lines represent the results obtained with the ap- 
proximate formula and ABACUS, respectively. The time scale is 

r = NTrh^/{2mgf). 

the total energy. The elementary excitations are quasiparticle- 
quasihole excitations in the Bethe-ansatz wave function ll23ll 
and may have mixed phonon and soliton character. Assum- 
ing that these have little effect on successive transitions, we 
estimate the decay of the center-of-mass velocity v by the 
classical equation Nmv ~ — Fv(i'), where is given by 
Eqs. ([1]) and (|4]i. This was integrated numerically and the 
result is shown in Fig. H) At the initial supersonic velocity, 
where the drag force is saturated (see Fig. [T]) the supercur- 
rent experiences constant deceleration. For v < c the drag 
force decreases and the deceleration slows down. For the 
TG gas we find an analytical solution for exponential decay 
v{t) ~ woexp(— t/r) for vq ^ vp. In the weakly-interacting 
regime, the decay becomes slow compared to experimental 
time scales. 



V. CONCLUSION 

Concluding, although the ID Bose gas with finite repulsive 
inter-particle interaction shows superfluid phenomena of the 
equilibrium type, we show that in general its ability to sup- 
port dynamic superfluid phenomena such as persistent ring 
currents is limited to a regime of very weak interactions; for a 
periodic potential, braking the gas, the persistent currents can 
be observed even in the TG regime at specific values of the 
velocity and density. 
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